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SIMULATION OF A CONTROLLED AIRFOIL WITH JETS 


BRIAN G. ALLAN *, MAURICE HOLT* , AND ANDREW PACKARD* 

Abstract. Numerical simulations of a two-dimensional airfoil, controlled by an applied moment in pitch 
and an airfoil controlled by jets, were investigated. These simulations couple the Reynolds- averaged Navier- 
Stokes equations and Euler’s equations of rigid body motion, with an active control system. Controllers 
for both systems were designed to track altitude commands and were evaluated by simulating a closed- 
loop altitude step response using the coupled system. The airfoil controlled by a pitching moment used an 
optimal state feedback controller. A closed- loop simulation, of the airfoil with an applied moment, showed 
that the trajectories compared very well with quasi-steady aerodynamic theory, providing a measure of 
validation. The airfoil with jets used a controller designed by robust control methods. A linear plant model 
for this system was identified using open- loop data generated by the nonlinear coupled system. A closed- 
loop simulation of the airfoil with jets, showed good tracking of an altitude command. This simulation also 
showed oscillations in the control input as a result of dynamics not accounted for in the control design. This 
research work demonstrates how computational fluid dynamics, coupled with rigid body dynamics, and a 
control law can be used to prototype control systems in problematic nonlinear flight regimes. 

Key words. Navier- Stokes equations, numerical simulation, feedback control, nonlinear flight regime, 
jets, coupled fluid dynamics, rigid body dynamics, and controls 

Subject classification. Applied Numerical Mathematics 

1. Introduction. In nonlinear flight regimes, the interaction between fluids, body dynamics, and con- 
trols can critically effect the performance of an aircraft. By coupling these disciplines one can computationally 
investigate the problematic nonlinear portions of the flight envelope. However, simulation models which rely 
on empiricism or linearity assumptions may give misleading results when applied to nonlinear flight regimes. 
To capture the nonlinear nature of the flow field, higher-ordered models need to be considered. In this 
investigation, the Reynolds- averaged Navier-Stokes (RANS) equations are coupled with Euler’s equations of 
motion, along with an active control system. 

The application which will provide the impetus for this work is the controlled landing of a high- 
performance powered- lift aircraft. This flight regime is currently problematic due to the difficulties in 
obtaining an adequate representation of the stability derivatives. In turn, use of a poor plant description to 
design the control system can lead to the loss of aircraft and pilot. The results of this research would begin 
to provide a control system designer with a means of computationally testing the control system in these 
problematic nonlinear flight regimes, thus helping to avoid these costly errors. 

Previous efforts have validated the component problems of unsteady fluid dynamics [1], the specified 
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Fig. 1. Overall coupled system approach. 


trajectory of a descending powered- lift aircraft [2], and the coupled fluids/body dynamics problem [3]. The 
coupled fluid, dynamics, and controls problem has been previously investigated for a three-dimensional store 
separation problem [4]. The case demonstrated here uses a similar methodology to describe the fluid and 
body dynamics by using a diagonalized implementation of the RANS equations in an overset mesh framework 
and Euler’s equations of motion to describe the rigid body response to aerodynamic loads. This study is 
different from the store separation problem as the interaction between fluid dynamics, rigid body motion, 
and controls, for the store problem was small. This was a result of the large ejection force on the missile, 
which dominated the trajectory of the store. Here, the interaction between fluids, rigid body dynamics, and 
controls is much more complex. The problem explored here will demonstrate the potential of coupling these 
disciplines for prototyping control systems in nonlinear flight regimes. 

2. Approach. The coupling used to solve this problem is shown schematically in Fig. 1 . This simulation 
starts with a converged fixed grid solution. The integrated aerodynamic loads and body states are passed 
from the flow solver to the six degrees-of- freedom, rigid body dynamics code. The new body position is then 
integrated one time step given the mass, moment-of-inertia, and applied forces. Any kinematic constraints 
on the body are also applied to the body motion at this time. The new body state is then passed to the 
controls module which generates a control input given the commanded altitude. Finally, the overlapped 
airfoil and jet grids are repositioned on the background mesh. This repositioning of the grids requires that 
the intergrid communication between the overlapped grids be re-established. This entire process is repeated 
every time step until the simulation is complete. 

2.1. Flow Solver. The flow field is computed by solving the RANS equations, using the diagonal 
scheme of Pulliam and Chaussee [5] in the grid framework of Steger et al [6]. The equations are inte- 
grated through Euler implicit time marching and second-order spatial differencing with viscous wall con- 
ditions specified as no-slip, zero normal pressure gradient, and adiabatic. Information transfered between 
overset mesh boundaries was implemented using trilinear interpolation of the dependent variable vector, 
Q = [p,pu,pu, pw, e ] T . The flow solver cost is 7ps/cell/step using a single processor of a Cray C90. Turbu- 
lence was modeled using the algebraic model of Baldwin and Lomax [7], as implemented by Renze [8]. 

2.2. Domain Decomposition. Computation of the loads generated by a dynamic vehicle requires 
accurate representation of the geometry and flow field. For the cases depicted in Fig. 2 and 3, the domain 
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was discretized using a curvilinear airfoil and nozzle grids overset on a stretched Cartesian background 
mesh. Overlapped grid topologies of this type allow representation of complex geometries and movement of 
bodies without regeneration of component grids [2, 4, 9]. Since grid boundaries are not required to match 
neighboring grids, robust algebraic [10] and hyperbolic [11] grid generators can be used to obtain spacing 
and orthogonality control. 

The exchange of flow information between the zones of an overlapping grid system are computed using 
a domain connectivity function. The donor- receiver relationship between the grids was established at each 
time step using an efficient search technique [12]. The relative cost of the intergrid communication to flow 
solver expense is dependent on the ratio of intergrid boundary points to volume points. 



Fig. 2. A closeup view of the airfoil grid overset on a Cartesian background grid. 

The overset mesh topology used for the airfoil with jets is shown in Fig. 2 and 3. Here, the intergrid 
boundaries are comprised of the outer boundary of the airfoil and nozzle grids, as well as the hole boundary 
in the background and airfoil grids. Along the hole boundary in the background grid, Q is trilinearly 
interpolated from the solution field of the airfoil grid. Similarly, the outer boundary points of the airfoil grid 
are interpolated from the dependent variable vector of the background grid. This exchange of information 
establishes the intergrid communication between the grid zones for steady or unsteady fixed grid problems. 
Note that the boundary surface points of each implicit grid zone must be specified either by boundary 
conditions or through communication with other zones. For problems in which zones move relative to 
each other, as is the case for aero dynamically controlled vehicles, the intergrid communication must be 
re-established at each time step. 

Note that Fig. 3 shows the bottom nozzle grid at the trailing edge of the airfoil with an outline of the top 
nozzle grid. The upper part of the bottom nozzle grid, which extends above the airfoil surface, is fictitious 
and does not interact with the flow field. These nozzle grids are used to allow the flow to develop as the 
mass flow rate is changed by the active controller at the inflow boundary. 

2.3. Rigid Body Dynamics. Rigid body motion is described by a conservation of linear and angular 
momentum using coordinates fixed to the center of gravity of the body and aligned with its principal axes [3] . 
Euler parameters are used to describe the rotation of the body from its initial position. These parameters are 
integrated according to the rotational body dynamics, updated, and stored for each grid [12, 13]. Kinematic 
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Fig. 3. 


The bottom nozzle grid at the trailing edge of the airfoil with an outline of the top nozzle grid. 


constraints can be imposed during the trajectory for restricted degrees-of-freedom simulations. In addition, 
the assumption of rigid body dynamics eliminates the need to store the component grids for all time, 
inasmuch as the Euler parameters may be used to compute grid attitude from the initial position. Note that 
the time scales of fluid motion are typically many orders smaller than the rigid body motion, making this 
loose coupling between the body motion and flow solver adequate. 

2.4. Control Law. The airfoil controlled by an applied pitching moment used an optimal state feedback 
control law. Integral control was used to improve tracking performance of altitude commands. The altitude 
of the airfoil was controlled by adding an applied moment to the aerodynamic forces, changing the attitude 
of the airfoil. This change in attitude produces aerodynamic lift, accelerating the airfoil vertically. The 
controller was designed using a linear time invariant system, based on quasi-steady aerodynamic theory. 
The feedback gains were then calculated using a Linear Quadratic Regulator (LQR) method, resulting in 
optimal feedback gains for the linearized model. 

In the second problem, the applied moment was replaced by two reaction control jets near the trailing 
edge. Due to the nonlinear nature of this system, a robust control design approach was used [14]. The 
linear plant model used in the control design was developed using input/output data generated by the 
coupled system. By using robust control design methods, a nominal linear model of the airfoil system can be 
parameterized by a measure of uncertainty. This parameterized uncertainty introduces unmodeled dynamics 
not accounted for by the nominal model. Therefore, by increasing the uncertainty in the model, the control 
design becomes increasingly robust to unmodeled dynamics not present in the linear nominal model. 

3. Results. The coupled system methodology described above was applied to a two-dimensional airfoil 
with and without jets. Both systems had controllers designed to track altitude commands. The closed-loop 
response of the airfoil without jets, using the coupled system, was compared to a quasi-steady aerodynamic 
model. The tracking performance of the airfoil with jets was evaluated from a closed-loop simulation of the 
coupled system. 

3.1. Airfoil Controlled by Applied Moment. This system consists of a NACA 64A010 airfoil with 
an applied pitching moment added to the aerodynamic forces. The airfoil has two degrees- of- freedom, one in 
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Fig. 4. Aerodynamic forces and notation . 

the vertical direction and the other in the pitch direction. A comparison of the closed- loop response of the 
coupled system was made by replacing the RANS equations with a quasi-steady aerodynamic model. This 
comparison will provide some measure of validation for the coupled system. 

3.1.1. Aerodynamic Model. The aerodynamic lift and moment coefficients for the airfoil are modeled 
using a quasi-steady formulation [15]. These coefficients are expressed as 

CUt) = C Lo + C Lo a(t ) + C U ~a(t) + C Lq l q{t ) 

Cm(U) = Cm„ + CM a a(t) + Cjji y &(t) + Cm 9 yQ(t) 

where Cl and Cm are the coefficients of lift and pitching moment about the quarter chord location respec- 
tively. The aerodynamic lift and moment in terms of Cl and Cm were 

m = C L {t) Q/>F 2 s) , M{t) = C M (t ) QpV 2 Sc) 

The stability derivatives Cm^ and Cm 9 were evaluated from steady state solutions of the RANS 

equations. The Cl & and Cm a stability derivatives were estimated from unsteady flow solutions of a plunging 
airfoil at a reduced frequency of k = 0.7. The terms, Cl 0 and Cm 0 represent the lift and moment coefficients 
evaluated at a = a = q = 0. All stability derivatives were evaluated for a free stream Mach number of 0.3, 
Reynolds number of 32 * 10 6 , and a chord length of 4.6 m. The values calculated for the quasi- steady model 
are given in Table 1 . 


Table 1 

Quasi-Steady Coefficients 


Cl o 

0.00 

C Mq 

0.00 

Cl 0 

6.28 

Cm, 

0.00 

CL q 

1.44 

C Mq 

-0.36 

Cl, 

0.78 

Cm, 

-0.35 


3.1.2. Linearized Dynamic Model. The equation of motion used for the vertical direction, mz = 
F z , can be expressed as 

(1) mz(t) = L(t) cos (#(£) — a(t)) 

where gravitational and drag forces are neglected. Since the lift, L(t) is perpendicular to the total velocity 
of the airfoil, as shown in Fig. 4, it was multiplied by cos(0 — a) to obtain its vertical component. 
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Fig. 5. A block diagram of the closed-loop system. 

The equation of motion in the pitch direction, I yy 0 = YL Af y , can be expressed as 

(2) lyyQit) = Af (t) + L(t) cos (a(£)) ^ + M u (t) 

where M u was the applied pitching moment produced by the controller. 

The quasi-steady aerodynamic lift and moment equations were substituted into Eq. (1) and (2) and 
linearized about the zero state. These equations were expressed as the linear time invariant system 

(3) xi (t) = Ax xi (t) 4- 
where xi = [0, 6 , z, z] T . Numerical values for A\ and B\ were 


Ax = 


These matrices were calculated using m — 62.8 kg , I yy = 386 kg • m 2 , — 0.3, and the coefficients in 

Table 1. 

3.1.3. State Feedback Design. The state feedback controller used in this study, incorporates integral 
control to improve tracking performance. Figure 5 shows a block diagram of the state feedback controller 
used. The objective of this control design is to find the optimal state feedback gains such that the altitude, 
z(t), tracks a commanded altitude path, z c (t ). Therefore, the error was defined as 

e(t) = z(t) - z c (t) - Dxi(t) - z c (t) 

Using state and integral feedback, the input u was 

u(t) = GiXi(t) G 2 J G(t)dt 

where G \ and G<i are the constant feedback gains. Now let, ± 2 (£) = e(t), where 
(4) ± 2 (t) = Dxi(t) - z c {t) 

Combining the hnear systems from Eqs. (3) and (4) results in 


Xi 

x 2 


Ax 0 
D 0 


xi 

X2 


-f 


Bx 

0 


u + 


0 

-I 
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This can be expressed as one linear system 


x(t) = Ax(t) + Bu(t) + Fz c (t) 


where x(t) = [xi(£) X 2 (t)] T . The input u(t) for state feedback becomes 


u{t) = 


' G x ' 

x i(t) 

g 2 

_ X2 (<) 


= Gx(t) 


Optimal feedback gains were calculated by finding the input, u(f), which minimizes the quadratic integral 
cost function 


J (tx(t)) 



(x T (t)Qx(£) + u T (t)Ru(t)) dt 


where Q and R are positive-semidefinite symmetric matrices. The quantity, x T (f)Qx(t), is a measure of how 
far the states deviate from the zero state at time t. The matrix Q is a diagonal weighting matrix, where 
each diagonal element corresponds to a penalty on one of the states. The term, u T (t)Ru(t ), accounts for the 
amplitude of the input used to bring the system to the zero state. This optimal control design is known as 
the LQR problem. 

The solution to this optimization problem is well known [16] and it can be shown that the input which 
minimizes J(u), given the system x(t) = Ax(£) + Bu(t) and the initial condition x(0) = 0, is 


u opt (t) = -R~ 1 B T Xx(t) 


Here, X is the solution to the algebraic Riccati equation 

A r X + XA- XBR~ l B T X + <2 = 0 

where A — BR~ l B T X is stable. Thus, the optimal state feedback gains are 

G = —R~ 1 B t X 


The optimal feedback gains G were calculated by first choosing initial values for the weights Q and R. 
After calculating the feedback gains G , the closed- loop response was simulated by solving Eqs. (1) and (2) 
coupled with the quasi-steady aerodynamic model. These coupled equations were integrated in time using 
a fourth order Runge-Kutta scheme. To assess the performance of the controller, a step in the commanded 
altitude of one chord length was simulated. The desired transient response for this simulation was a settling 
time of 1.0 s with minimum overshoot. If this performance criteria was not achieved, the weights were 
adjusted to increase or decrease the penalty of the individual states. This change in the weights resulted in a 
optimal controller with a new transient response. This cycle was iterated upon until the desired performance 
was achieved. The final iteration produced the feedback gains 

G = [—193 - 13.1 - 2.43 1.07 - 3.00] - 10 3 

3.1.4. Airfoil with Applied Moment Simulation. A step command in altitude of one chord length 
from the starting position was simulated. Figure 6 shows states of the controlled simulation for the linear 
aerodynamic model and the nonlinear coupled system response. These two simulations match up very well 
as would be expected in this linear flight regime. This comparison provides a measure of validation for the 
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Fig. 7. Instantaneous Mach contours during the ascent of the airfoil to the commanded altitude. 


nonlinear coupled system. Instantaneous Mach contours of the simulation for the ascent of the airfoil to the 
reference altitude, is shown in Fig. 7 at four different times. 

In this simulation, the applied moment pitches the airfoil up producing aerodynamic lift. This lift 
accelerates the airfoil in the vertical direction, moving the airfoil to the commanded altitude. As the airfoil 
starts to reach the commanded altitude, the controller produces a downward pitching moment. This applied 
moment pitches the airfoil down which reduces the lift generated by the airfoil, slowing down its ascent. 
Figure 6 shows a small amount of overshoot of the commanded altitude, as predicted in the control design 
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using the quasi-steady model. This overshoot was also seen in the RANS coupled simulation. A slight 
difference in the control moment generated can be seen in Fig. 6 where the quasi-steady model predicts a 
smaller control input. 



Fig. 8. Two-dimensional airfoil with jets problem. 


3.2. Airfoil Controlled by Jets. In this problem, the applied moment is replaced by two reaction 
control jets near the trailing edge of the airfoil. The center of gravity is placed at the quarter chord location 
as shown in Fig. 8, increasing the stability of the system. 

It should be noted that the problem of resolving the complex flow of a two-dimensional jet in a cross 
flow was not addressed here, since the purpose of this investigation was to demonstrate the analysis of a 
control system in a nonlinear flight regime. Therefore, an approximation to the jet flow was made by using 
coarse grids which do not resolve the complex flow associated with jets. However, it maybe argued that the 
approximation made here is satisfactory for control design and analysis purposes as long as the flow being 
resolved captures the general dynamics of the system. 

3.2.1. Plant Model Design. The highly nonlinear aerodynamic effects, from the jets at the trailing 
edge of the airfoil, make it very difficult to model. Figure 9 shows the aerodynamic lift and moment time 
histories for a fixed airfoil with a constant mass flow rate from the bottom jet. The aerodynamic forces have 
a large oscillatory behavior produced by the unsteady flow generated behind the jet. The large magnitude 
of these oscillations were due to the jet being two-dimensional. The power spectral density of the Cl time 
history from tU^/c = 0 to 10 is shown in Fig. 10. This figure shows that the lift generated by a constant mass 
flow rate to the bottom jet does not have a dominate frequency. The approach used to model the dynamic 
response of the airfoil was to perform a system identification using input/output data. This input/output 
data was generated from two open-loop simulations using the coupled system. Note that a turbulence model 
was not used for this problem as it was not valid for this complicated unsteady flow field. 

The relationship between the input u and nondimensional mass flow rate to the top and bottom jets is 
shown in Fig. 11. A positive input u would correspond to a mass flow rate to the bottom jet and a negative 
u for the top jet, where u has a magnitude limit of 0.3. 

The variables Ur and u B are functions of the input u and are defined in the following way: 



[° 

u<0 


r 0.3 

u < -0.3 

Ut = < 

1 u 

0 < a < 0.3 , 

u B = < 

- 

-0.3 < u < 0 


[o,3 

0.3 < u 


lo 

0 < u 


A delay for the mass flow rate command and the actual mass flow rate applied to the nozzles is modeled 
by 1 /(t s + 1), as shown in Fig. 11, where r = 4.5 * 10“ 4 . 

The nondimensional mass flow rates, ra\ ov and were applied to the top and bottom nozzle grids 
respectively. These mass flow rates were used to change the inflow boundary conditions on the nozzle grids 
shown in Fig. 3. Note that a value of m* =0.1 would corresponds to a mass flow rate ofm= 13.7 kg/s at 
the inflow boundary. 
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Fig. 9. Time history of Cl and Cm for a fixed airfoil with the bottom jet turned on for a nondimensional mass flow rate 
of 0.3 from tUoo/c =0 to 10. Time is nondimensionalized in terms of chord distances flown. 



0 20 40 60 80 100 120 140 

Frequency, Hz 


Fig. 10. The power spectral density of the Cl time history shown in Fig. 9 from tUoo/c — 0 to 10, where Uoo — 102 m/s 
and c — 4.6 rn. 

The state space model developed has an input u, with outputs i and 9. The altitude z was found by 
integrating i. The pitch rate was produced by restructuring the linear system to obtain 6. This was done 
by first looking at the output equation 



Thus the pitch rate could be expressed as 

9 = c 2 x = c 2 (Ax 4- Bu) 

Restructuring the output equation to include 9 resulted in 



This system, P{s) y is shown in Fig. 12 where the altitude was found by integrating the vertical velocity. 

The input u and states x, for the open-loop simulations, are shown in Fig. 13 where u(£) was a series 
of step commands. The first simulation is shown in Fig. 13a and the second in Fig. 13b. The model was 
identified using both sets of input /output data simultaneously. A comparison of the coupled simulation 
and linear model is shown in Fig. 13. This comparison shows how the simple two state model captures 
the general behavior of the coupled system for the given inputs. The model compares relatively well for 
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m t0 p(s) 

m'boi(s) 


u(s) 


R(s) 

U T (S) 

1 



xs + 1 


1 

Ub(s) 

TS + 1 



Fig. 11. Block diagram showing the relation between the input u and the nondimensional mass flow rates to the top and 
bottom jets in the Laplacian domain. 


Z(s) 

Z(s) 

0(s) 

0 ( 8 ) 

Fig. 12. State space model with integrator and added 6 output. 

the second simulation and for the beginning of the first simulation. In the first simulation, the model and 
coupled response start to diverge near t = 0.1 s. This could be a result of the larger attitude and pitch rate 
as compared to the second simulation. 

3.2.2. Robust Control Design. The model used for the robust control design consists of the nominal 
model P(s) and the elements, 7 and A, which parameterize the uncertainty in the model. This system is 
shown in Fig. 14 where 7 was a scalar and determines the amount of uncertainty in the model. The transfer 
function A is assumed stable and unknown, except for the norm condition, || A||oo < 1. The interconnection 
for the inputs w and outputs v to the A block are shown in Fig. 14a. The perturbed system with A connected 
is shown in Fig. 14b. 

The block diagram in Fig. 15 shows the performance weights and disturbances used for the robust control 
design [14]. The W m0( iei block was the ideal model response, which was represented by a well damped second- 
ordered system. The performance weight W cm( i defines the magnitude and frequency of the input commands 
to be tracked. W actx and W act2 were penalty weights on the input u and its rate u. The difference in the 
altitude and the ideal altitude response was weighted by W per f X . The attitude, pitch rate, and vertical 
velocity were weighted by W per f 2 . Noise was added to the system by inputs cfo to d 5 and shaped by W snois 
which sets the magnitude and frequency range. 

The performance weights and plant G in Fig. 15 were reduced to the block M shown in Fig. 16. The 
controller K was designed so that the perturbed system was stable and the transfer function from d to e 
satisfies 

||T ed (M,^,A)|| oo <l 

for all stable perturbations, A satisfying ||A||oo < 1 

The controller K was determined by using ji synthesis, which was approximated by an iterative method 
known as U D — K iterations” [14]. The original controller K, which had 68 states, was reduced to 17 states 
by performing a balanced realization [17]. This reduced controller was then used for the coupled RANS 
simulation. 

3.2.3. Airfoil with Jets Simulation. The input y to the controller K was a vector containing the 
commanded altitude and the states, where y = [z c ,z, z, 6 , 0\. The controller takes these inputs and generates 
the signal u(t), which was passed through the block shown in Fig. 11. This produces mass flow rates for the 
inflow boundary on the nozzle grids. 
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Fig. 13. Two open-loop simulations of the nonlinear coupled system with the linear input/output model response. 





Fig. 14. Uncertainty model for P(s). a) Uncertainty block A not connected where 7 determines the amount of uncertainty; 
b) Plant G(s) with full uncertainty block A. 


Performance of the controller was evaluated by analyzing the closed-loop response of the airfoil to a step 
input in the commanded altitude. Figure 17 shows the states for the simulations using the the linear and 
nonlinear models. The altitude trajectory for the linear model was very good and was nearly identical to 
the ideal model response. The altitude trajectory for the coupled system was also good. It showed that the 
controller was able to track the desired altitude trajectory within the prescribed tracking performance. 

The nonlinear simulation also revealed oscillations in the control jets not shown in the linear model 
simulation. The difference between these two simulations show how the nonlinear dynamics were not captured 
by the linear plant model. These oscillations are induced by the controller, which is trying to force the 
nonlinear airfoil system follow the desired altitude response of a linear system model, W^odeh The frequency 
of these oscillations are thought to be partly related to the altitude tracking performance, W perf 1 , used to 
design the controller. By increasing the tracking performance, the controller will try to keep the airfoil even 
closer to the desired altitude path, increasing the frequency of the oscillations. 

The instantaneous Mach contours for the coupled simulation are shown in Fig. 18 at four different times. 
This figure reveals the complicated flow pattern generated by the top and bottom jets during the closed-loop 
simulation. The snapshot at t = 0.02 s was taken at the beginning of the simulation where the top jet was 
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Fig. 15. Block diagram of performance weights for robust control design. 
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Fig. 16. Closed-loop system with performance weights and plant model G contained in M . 


turned on by the controller. This top jet produces aerodynamic forces pitching the airfoil up. As the airfoil 
pitches up, aerodynamic lift is generated and begins to move the airfoil vertically. The instantaneous Mach 
contours at t = 0.4 s show the unsteady wake behind the airfoil as a result of the oscillations in the top and 
bottom jets. A closeup view of the trailing edge can be seen in the last two snapshots of Fig. 18. This view 
shows the unsteady flow behind the top jet where a low pressure region has developed. At t = 1.2 s the 
airfoil has reached the commanded altitude. 

4. Conclusions. The analysis of an aircraft controls system for nonlinear flight regimes is complicated 
by the interaction of fluid dynamics, rigid body dynamics, and the control system itself. Of particular interest 
is the landing of a powered-lift aircraft which has a strong interaction between the lifting jets and ground 
plane. Towards this goal of analyzing this system, the coupling of the Reynolds- averaged Navier-Stokes 
equations and Euler’s equations of rigid body motion with an active control system, was performed. As a 
demonstration problem, the altitude control of an airfoil with and without jets was investigated. 

The airfoil without jets was controlled by an applied moment and used an optimal state feedback 
controller with integral control. The control design used a Linear Quadratic Regulator method with a 
linear time invariant system based on quasi- steady aerodynamics. The stability derivatives for this model 
were computed using steady and unsteady solutions to the RANS equations. A step in the commanded 
altitude was simulated by the RANS coupled system which compared well with the quasi- steady model. 
This comparison provided a measure of validation for the nonlinear coupled system. 

The airfoil controlled by jets used a robust control design method. The linear plant model, used in 
the control design, was developed by performing a system identification on the input/output data. This 
input/ output data was generated by two open-loop simulations of the coupled system. Simulation of the 
closed-loop response of the airfoil with the coupled system, demonstrated that the airfoil was able to track 
the desired altitude trajectory within the specified performance. It also showed oscillations in the control 
jets that were not predicted by the linear model simulation. 
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Fig. 18. Instantaneous Mach contours during the ascent of the airfoil to the commanded altitude. 


REFERENCES 

[1] C. A. Atwood and W. R. Van Dalsem, Flowfield simulation about the SOFIA airborne observatory , 

AIAA J. of Aircraft, 30 (1993), pp. 719-727. 

[2] K. Chawla and W. R. Van Dalsem, Numerical simulation of STOL operations using thrust- 

vectoring, AIAA Paper 92-4254, 1992. 

[3] R. L. Meakin and N. Suhs, Unsteady aerodynamic simulation of multiple bodies in relative motion , 

AIAA Paper 89-1996, 1989. 

[4] C. A. Atwood, Computation of a controlled store separation from a cavity, AIAA J. of Aircraft, 32 

(1995), pp. 846-852. 

[5] T. H Pulliam and D. S. Chaussee, A diagonal form of an implicit approximate-factorization algo- 

rithm, J. Comput. Phys. 39 (1981), pp. 347-363. 

[6] J. L. Steger, F. C. Dougherty, and J. A. Benek, A Chimera Grid Scheme , in Advances in Grid 


14 









Generation, K. N. Ghia and U. Ghia, eds., American Society of Mechanical Engineers, New York, 
1993, pp. 59-69. 

[7] B. S. Baldwin and H. Lomax, Thin-layer approximation and algebraic model for separated turbulent 

flows , AIAA Paper 78-257, 1978. 

[8] K. Renze, P. Buning, AND R. Rajagoplan, A comparative study of turbulence models for overset 

grids , AIAA Paper 92-0437, 1992. 

[9] R. L MEAKIN, Moving body overset grid methods for complete aircraft tiltrotor simulations , AIAA Paper 

93-3350, 1993. 

[10] J. P. Steinbrenner, J. R. Chawner, and C. L. Fouts, A structured approach to interactive 

multiple block grid generation , in Advisory Group for Aerospace Research & Development, Confer- 
ence Proceedings No. 464, Applications of Mesh Generation to complex 3-D Configurations, 1989, 

pp. 8.1-8.12. 

[11] W. M. Chan and J. L. Steger, Enhancements of a three-dimensional hyperbolic grid generation 

scheme , Appl. Math, and Comput. 51 (1992), pp. 181-205. 

[12] R. L Meakin, A new method for establishing inter-grid communication among systems of overset grids , 

AIAA Paper 91-1586, 1991. 

[13] R. L Meakin, Computations of the unsteady flow about a generic wing /pylon/ finned- store configuration , 

AIAA Paper 92-4568, 1992. 

[14] A. Packard, J. Doyle, and G. Balas, Linear, multivariable robust control with a p perspective , 

ASME J. of Dynamic Systems, Measurement, and Control 155 (1993), pp. 426-437. 

[15] M. Tobak and L. B. Schiff, Aerodynamic mathematical modeling - basic concepts , in Advisory Group 

for Aerospace Research & Development, Lecture Series No. 114, Dynamics Stability Parameters, 
1981, pp. 1.1-1.32. 

[16] H. Kwakernaak and R. Sivan, Linear Optimal Control Systems , John Wiley & Sons, New York, 

1972. 

[17] K. Glover, All optimal hankel-norm approximations of linear multivariable systems and their L ^ 

error bounds , International J. of Control, 39 (1984), pp. 1115-1193. 


15 



Form Approved 
OMB No. 0704-0188 


REPORT DOCUMENTATION PAGE 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports 1215 JefFerson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget. Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

October 1997 


3. REPORT TYPE AND DATES COVERED 

Contractor Report 


4. TITLE AND SUBTITLE 

Simulation of a Controlled Airfoil with Jets 

5. FUNDING NUMBERS 

C NAS1-19480 
WU 505-90-52-01 

6. AUTHOR(S) 

Brian G. Allan 
Maurice Holt 
Andrew Packard 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Computer Applications in Science and Engineering 
Mail Stop 403, NASA Langley Research Center 
Hampton, VA 23681-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

ICASE Report No. 97-55 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-2199 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR- 20 1750 
ICASE Report No. 97-55 

11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Dennis M. Bushnell 
Final Report 

To appear in the Journal of Guidance, Control, and Dynamics 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category 64 
Distribution: Nonstandard 
Availability: NASA-CASI (301)621-0390 

12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

Numerical simulations of a two-dimensional airfoil, controlled by an applied moment in pitch and an airfoil controlled 
by jets, were investigated. These simulations couple the Reynolds-averaged Navier-Stokes equations and Euler s 
equations of rigid body motion, with an active control system. Controllers for both systems were designed to track 
altitude commands and were evaluated by simulating a closed-loop altitude step response using the coupled system. 
The airfoil controlled by a pitching moment used an optimal state feedback controller. A closed-loop simulation, of 
the airfoil with an applied moment, showed that the trajectories compared very well with quasi-steady aerodynamic 
theory, providing a measure of validation. The airfoil with jets used a controller designed by robust control methods. 
A linear plant model for this system was identified using open-loop data generated by the nonlinear coupled system. 
A closed- loop simulation of the airfoil with jets, showed good tracking of an altitude command. This simulation 
also showed oscillations in the control input as a result of dynamics not accounted for in the control design. This 
research work demonstrates how computational fluid dynamics, coupled with rigid body dynamics, and a control 
law can be used to prototype control systems in problematic nonlinear flight regimes. 


14. SUBJECT TERMS 

Navier-Stokes equations, numerical simulation, feedback control, nonlinear flight regime, 
jets, coupled fluid dynamics, rigid body dynamics, controls 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 




18. SECURITY CLASSIFICATIO 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATIO 
OF ABSTRACT 


15. NUMBER OF PAGES 

20 


16. PRICE CODE 

A03 


20. LIMITATION 
OF ABSTRACT 


tandard Form 298 (Rev. 

Prescribed by ANSI Std Z39-18 
298-102 













